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ABSTRACT 


The hig^ accuracy of contemporary observational techniques now requires 
a better physical and mathematical urderstandlng of the form of the tidal poten** 
tial acting on the satellite and a gradual reshaping of the computational methods 
currently in use. The dependence of elastic responses and tidal lags on the 
geographic position and on tidal frequency is introduced into the tidal potential 
on the surface of the Earth. By making use of Dirichlet's theorem we continue 
the tidal potential from the surface of the Earth into outer space. The final 
form of the disturbing tidal potential acting on the satellite is obtained in terms 
of the mean elements of the satellite and in terms of the standard arguments of 
the lunar theory. This form of the tidal potential is suitable for performing the 
analytical integration and for computation of the tidal effects over a long interval 
of time, because all agruments are linear or very nearly linear with respect to 
time. We interpreted the tidal effects in the motion of the satellite as perturba- 
tions in the orbital elements and as the small periodic variations in the co- 
efficients of the spherical harmonics in the expansion of the geopotential. In 
transforming the exterior tidal potential to a form containing the orbital elements 
it was convenient to resort to the expansion in terms of the generalized spherical, 
harmonics familiar from the theory of angular momentum in quantum mechanics. 
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BASIC NOTATIONS 


G - 
M - 
R - 


jo 


p', q', s' 


m 


t ^ 


a - 


't, F, D and F - 


0', - 
a', 8' - 


n' - 


cos 7' 


r 


» 


r 


M 


r 


«o 


m - 


gravitational constant, 
the mass of the Earth, 

the radius of the Earth, considered as a sphere, 

position vector of the Moon relative to the center of 
the Earth, 

the unit vector in the direction of r', 
the equatorial components of r'^ , 
the mass of the Moon, 

the mean distance of the Moon from the Earth, defined 
in such a manner that the constant part in the expansion 
of a'/r ' is equal to 1, 

the standard arguments of Brown’s lunar theory, 
the geogra^diic colatitude and the longitude of the Moon, 
the ri^ ascension and the declination of the Moon, 
the lunar tidal disturbing function, 
r® • r'o, 

position vector '• Sun relative to the center of the 
Earth, 

Ir1. 

the unit vector in the direction of r" , 
the equatorial components of r"® , 
the mass of the Sun, 
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a" > the mean distance of the Sun from the Earth, defined 

in such a manner that the constant part in the expansion 
of a"/r " is equal to 1, 


0* , - the geographic colatitude and longitude of the Sun, 

a", §" - the right ascension and the declination of the Sun, 
^ - the solar tidal disturbing function. 


0 - the sidereal time, 

0, \ - the geographic colatitude and the longitude of a point 
on the ^•'rth surface, 

r - position vector of the Satellite, relative to the center 
of the Earth, 


r° - the unit vector in the direction of r , 


a, b - the right ascension and declination of the Satellite, 


a. e, CO, Q, I 


the elliptic elements of the instantaneous orbit of the 
satellite. 


i - the true anomaly of the satellite. 


R - the unit vector normal to the instantaneous orbital 
plane of the satellite. 
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INTRODUCTION 


In the present work we Investigate the form of the exterior tidal potential 
acting on a satellite. The tidal potential can be defined as small 'periodic vari- 
ations in the geopotentiaL caused by the formation of the tidal buldge and by 
the redistribution of masses in the Earth's interior under the influence of luni- 
solar tidal attraction. In other words, the tidal effects in the motion of the 
satellite can be interpreted as the result of tidal periodic oscillations in the 
coefficients of the spherical harmonics in the expansion of the exterior geo- 
potential. Of considerable interest and importance are the long period oscil- 
lations caused by the regression of the lunar node along the ecliptic. The 
existence of such oscillations m the oblateness coefficient C 20 was first 
emphasized by Kozai (1965). At the present time the influence of the long period 
tidal effects on the coefficients in the expansion of geopotential lies on the limit 
of observability. However, with the further amelioration of the observational 
techniques and extension of observations over a longer interval of time we shall 
be obliged to take these effects into consideration. The variations in geopoten- 
tial produce in turn small tidal perturbations in the motion of the satellite. 

They can be clearly detected by the modern observational techniques (Newton, 
1968), (Smith et al., 1973). The high accuracy of contemporary techniques now 
requires a better physical and mathematical understanding of the phenomena and 
a gradual reshaping of the computational methods currently in use. The programs 
currently in use permit one to determine the average elastic properties of the 
Earth from the magnitude of the tidal effects in the motion of satellites. 
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By adjusting numerical values of parameters of the Earth's elastic response 
and of the tidal lags one can achieve a satisfactory representation of the satel- 
lites motlor It was found, however, that the numerical values of the mean 
elastic parameters, as determined from the motion of satellites, are different 
for different orbital inclinations (Smith et al., 1973), (Douglas et al., 1974). This 
clearly indicates the existence of lateral Inhomogeneltles and, possibly, asym- 
metry in the Earth's internal structure and of their Influence onthe tidal pertur - 
bations of a satellite. Until the problem of the influence of lateral Inhomogenei 
ties on the tidal oscillations is resolved, we are forced to use the information 
available on the Earth's surface, the parameters of the elastic response (Love 
coefficients) and tidal lags, and by making use of Dirichlet's theorem continue 
the tidal potential anal 3 rtically from the surface of the Earth into outer space. 

Thiis, for the present any further refinement in the theory of tidal per- 
turbations in the motion of satellites is tied, to considerable extent, to the 
measurements on a global scale of tides or tidal elastic responses on the sur- 
face of the Earth, directly or by remote sensing. In other words, the dependenct^ 
of elastic responses and tidal lags on the geographic position and on tidal fre- 
quency shall be introduced into the tidal potential. 

Original Love elastic parameters (Love, 1909) are the numerical coefficients 
attached to Legendre polynomials in the expansion of the tidal potential. Thus. 
Love formulation assumes that the elastic responses are the same for all tidal 


2 



constituents associated with a given Legendre polynomial, irrespectively of the 
tidal argument. Kaula (1969) and Balmino (1973) assumed that Love coefficients 
depend on geogriq;>hic position and are sufficiently smooth functions on the Earth 
surface so that they can be expanded into a series of spherical harmonics. 
Balmino assumed in addition a dependence of Love coefficients on the tidal fre- 
quency. This last assumption is in conformity with the results by Alterman et al. 
(1959). In the present exposition we follow Balmino and assume the dependence 
of Love numbers on the geographic position as well as on the tidal frequency. 

We found it convenient to use complex "Love numbers," because they carry the 
combined influence of elastic responses and lags on the motion of the satellite. 

It is necessary to remark that the introduction of Love coefficients pre- 
supposes the absence of non-linearities, i.e. the absence of noticeable (in the 
motion of a satellite) interactions between the tidal constituents (Kaula, 1969). 

Only the influence of long period tidal effects can be detected in the elements 
of artificial satellites. As a consequence the form of the exterior tidal potential 
can be greatly simplified, because all terms which involve the sidereal time or 
the mean anomaly of the satellite can be omitted from the expansion. This means 
also that the jppm^lei^ Love coefficients can be averaged along the parallels. The 
remaining dependence of the averaged Love numbers on latitude will produce 
the primary and the secondary tidal effects in the motion of the satellite. The 
secondary effects are of the same form and magnitude as some primary ones 
(Musen and Felsentreger, 1973). 
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Care must be taken against oversimplification of the theory and computa- 
tions by neglecting essential secondary effects or by confusing them with pri- 
mary ones. In transforming the tidal potential to a form containing the orbital 
elements it is convenient to resort to expansion in terms of spacial and surface 
generalized spherical harmonics. These harmonics are familiar from the theory 
of angular momentum in quantum-mechanics (Edmonds, 1960) (Rose, 1957), 
(Vilenkt‘ 1 : 1965). They also can be profitably applied in the theory of motion of 
celestial bodies and in theoretical seismology (Burridge, 1969), (Phinney and 
Bur ridge, 197^), (Smith, 1974). 

The F -coefficients in Kaula's (1969) elegant expansion of the tidal potential 
are intimately related to the generalized associated Legendre functions which are 
being used in the present work. 

Nnmnronr; recursive relations between the generalized spherical functions 
and their derivatives greatly facilitate the process of expansion of the tidal 
potential. Using the recursive formulas and an analytic programming language 
one can obtain a purely analytic expansion by computer or, by substituting the 
numerical value of the inclination (and of the eccentricity) at the very outset, 
one can obtain a semi- analytical expansion of the tidal perturbations, in the m 
of periodic series with purely numerical coefficients. 

The first option gives the solution in a most general form. The second one 
requires less computing time and makes the programming work easier. 
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We propose, as In our previous work (Miisen and Estes, 1972), to discard the 
expansion of the tidal potential in terms of the elements of the Moon, equatorial 
or ecliptical. In particular, the equatorial elements of the Moon are not linear with 
respect to time and therefore are not q\dte suitable for the integration in analytic 
cal or semi-analytical form, nor for the investigation of the tidal effects over a 
long interval of time. We propose instead a Doodson (1922) type expansion of 
the tidal potential in terms of the standard arguments of the lu 2 .ar theory 
with purely numerical coefficients. In Doodson type expansion all arguments 
are very nearly linear with respect to time. This facilitates integration of dif- 
ferential equations for the variation of eleme .ts and the search for the critical 
terms with low frequencies. In performing this expansion we can make use of 
the developments of the lunar rectangular coordinates directly available from 
Brown's lunar theory (with some modem corrections). All necessary programs 
to produce the semi-analytical expansions of the tidal potential and of the tidal 
perturbations in the elements of the satellites were prepared by R. Estes. 

(Musen and Ester, 1972). 

TIDAL DISTURBING POTENTIAL 


We start from the standard expansion on the surface of the Earth of the 



and £)olar 



tidal disturbing functions, respectively, with the Love elastic coefficients and 
tidal lags temporarily omitted. Introducing normalized surface harmonics 


Y. (4>, k) = 

w n 


y^I^PJ(co8*)exp(im\), 


(S) 


\.J 4 ; \) = (-irY;,(«. M 




with the normalization 


4t7 


k) dS = 


1 

2i + 


1 j ’ 


( 5 ) 


dS = sin i^di^dA., 


over the unit sphere, and making use of the addition theorem: 


P^(cos/)= r X) 


and with 


we obtain 


m 


- a' -e, 


n' = Gm 




exp(im0)Y^J|^ 


(6) 

( 7 ) 


( 8 ) 


where we set: 


Vm 



YJ (0'. a') 

L m 


( 9 ) 
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and a similar expression: 







( 10 ) 




m 



( 11 ) 


for the solar tidal disturbing function. The expressions (9) and (10) depend on 
the equatorial coordinates of the Moon and Sun, respectively. The lunar and 
solar tiieories currently in use give the expansions of coordinates, rectangular 
or polar, of both bodies in form of trigonometric series with the standard argu- 
ments ^ , F, D and r , and purely numerical coefficients. The accuracy of 
the expansions is more than sufficient for the computation of tidal effects in the 
motion of the satellite. 

The arguments 'i , Vt F, D and r are very nearly linear with respect to 
time and their use facilitates the integration of the differential equations for 
variation of elements and the long time prediction of tidal effects in the motion 
of a satellite (Musen and Estes, 1972). To facilitate the programming and the 
expansion into periodic series, we can express the "lunar" spherical functions 
Yo , a' ) in Maxwelliari form as polynomials in the equatorial components 

^ Hi 

P* , q\ s' of the unit vector directed toward the Moon from the center of the 
Earth. We have: 
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( 12 ) 





+ iq'r • 






(-1)" 


( 2 i - 2 k) ! j,-e-«- 2 k 

k ! ('t - k) ! (^ - m - 2 k)! 


k"0 

for m > 0. We can use (4) to obtain the expansion for m < 0. Expansions of 
p' , q' , s', a' A' and then of Ti can be obtained easily from Browns series 

'Vm 

(1919) of the lunar coordinates and parallax (with modem corrections) using 


existing computer programs (Musen and Estes, 1972). 

The resulting expansion of Tj is the sum of the expressions of the form 

Tj =N® exp(iwT). ( 13 ) 

't/m W 't-mV 


where 



V 


V=(jj» J2» 13' 34’ is) 


is a vector with integral components and 

T = (^, V, F. D. ! ) 


(10 


and are the numerical coefficients. It is convenient at this point to intro- 
duce Love coefficients and their dependence on geographic position and frequency. 
We replace (13) in the expansion of the tidal disturbing function by the expression 

(0. ^ ; T 'o , 

-VmV 

where 

Ko (0. \) = k. (p. X)exp[-i€. (<p, A)3. (15) 

■t-mW -t-mW 't'WV 

and ko . t. are the standard Love coefficients and tidal lags, Respectively. 

'tmV 'tmV 
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The lunar tidal disturbing function on the surface of the Earth is the sum 


of expressions of the form: 


"Omv 



exp(ime) \) 


(16) 


For the sake of compactness we shall temporarily suj^ress the dependence on 
V and write simply: 


R 


= Gm* — exp(im^) ^ **^1, 


Substituting the expansion 


+ ® cj ■ + p 


pTJj q*-p 


viiiere 


K. r J K. Y p Y* dS, 

'vmpq 4 t7 


(17) 


(18) 


(19) 


into (17), we obtain: 


exp(imtJ) T| 


♦ ® q ■ » p 

■ E Z! 

p • 0 q ■ ” p 


Y . 

'^mpq P*! 


( 20 ) 


Making use of (4) we have: 


K. =(-!)*» f K. Yp Y„ . dS. 

'I'inpq 4 tT 


( 21 ) 


Below we make use of theorems and formulas familiar from the theory of the 
representation of the rotation group and from the theory of angular momentum 
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in quantum-mechanics. We found it convenient in our work to follow Vilenkin's 
(1965) notation and symbolism as close as possible. 

Let x» 4> ^ Euler's angles defining a rotation g. The elements of 

the representation of order of g(the generalized spherical harmonics) have the 
form: 

'/') = exp[-i(mx + nv//)] P^^(cos 0), (22) 


m 


, n = - -i, 




where (cos0) are the geneialized associated Legendre functions. We define 
the angles x and in such a manner that (cos0) are always real. We have 
the relation: 

0, -) (23) 

with 


X; = X + 7^. 

From (23), (4) and making use of the Clebsch-Grordsn series: 


0 . </') = 


(24) 


(25) 




C(i^, s; jj, j 2 )C(^j, ^ 2 , s; kj, k 2 ) j j+i 2 < 0. 0). 


where the symbol C designates Clebsch-Gordan coefficients (Rose, 1957), we 
derive: 


Y ^ (- 1)' 

p,-q ' 


,m-q 


C('t, p, s; -m, +q)C('t, p, s; 0, 0) Y* 

9 y q n 


(26) 
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Substituting (26) into (21) we obtain: 


h 


mpq 


i* I'tijp I 

(-1)“ ^ C(iy p, s; -m, +q) C(^, p, s; 0, 0) 

s-PT-p! 


2p+ 1 ^ 

2s + 1 


( 27 ) 


where symbols .designate the coefficients in the e3q>ansion of the Love 

coefficient k . into a series of spherical harmonics, 


'vm s j 


2S + 1 r 

477 ^ 


Ko Y'.dS, 

't/m ■ J 


(28) 


In particular we have; 



viiere 



(29) 


(30) 


(31) 


is the value of k ^ averaged along the parallel. The analytical continuation of 

'vm 

the tidal potential (20) into extraterrestrial space is a harmonic function and 
therefore has the form; 

+ <D q ■ + p 
p • 0 q ■ ' o 

where is the lunar parallactic factor, 



R\p+1 


Yp^(0, \). (32) 
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In application to the motion of a satellite we set in (32): 


X. = a - 0, 



8 , 


and we obtain 


Q 


'tm 


_ Gm' + 


R 






( 33 ) 


+® Q * P 

Z] Z] (t)' ' “>■ 

p = 0 q* -p 

where j6 is the satellite's parallactic factor. 



We are not interested in the tidal effects with period of one day or less. There- 
fore only the terms for which m = q are of any practical importance and shall 
be kept in (33). All remaining terms can be discarded and we can set: 



where 


o; - 

<-mp R 


_ Gm' ,'t+i 


j3' 


mpm 


'Vm 





a), 


(35) 


and from (27) we have: 

K+p| 

s» [T^pi 

and is given by the combination of (30) and (31). 
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The form 







pm 


( 0 . 


of ( 35 ) indicates that 


8C, 

ra 


m/ 

M 



h 


mpm 


can be interpreted as perturbations in the coefficient 
harmonic 



'pm 


exp(im0) 


of the spherical 





in the expansion of the exterior geopotential. 

With this interpretation the coefficients in the geopotential become chang- 
ing periodic fimctions of time. For perturbations in the coefficient of the zonal 
harmonic (R/r)'*^^ Yp^ we have; 


'IpO M 




In particular for the perturbations in the oblateness -coefficient C 20 we have: 




An estimate of the pertiirbations in C20 is given by the formula: 


8C„ = i/3'^k, » P,(s') 


M 


»\3 
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where = 0.3, the approximate value of the Love number, and (a'/r')^ P 2 (s') 
shall be expanded into a periodic series with the arguments of the lunar theory. 
The terms in which contain only the longituie of the lunar ascending node 
- F + D + r produce "very long period" effects with the periods of approxi- 
mately 18, 9, , . . years. 

At the present time the influence of such '*very long period" terms, with the 
periods of the order of years, on the coefficients in the expansion of the geopoten- 
tial (or on the motion of satellites) lies on the limit of observability, but with the 
further improvement of the observational techniques and extension of observa- 
tions over a longer interval of time the importance of these terms might increase 
and we will be obliged to take them into consideration, especially in C 20 • 
Application of the addition theorem 

k = 

tl(g2> (37) 

k-^e 

for the generalized spherical harmonics, as defined by the rotations g^ and 
respectively, probably represents the shortest way to obtain the expansion of 
the tidal disturbing function (35) in terms of the angular elements of the satellite's 
orbit. The transition from the frame as defined by the unit vectors (r , R x r® ,R) 
to the system associated with the equator involves the rotations: 

gr 0=1, 0 = =f+w+— , 
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g^: X = '^« ‘^=■5- 


'7 


8jR2= X = <P--2~^^ 


From (37) and taking (22) -(23) into account we obtain: 

^ Ko{^’ i ' -) ■ “p{‘ 


k--p 




From (22) and making use of relations 


PEk,oW = 


''Eo<’‘> V(T7§T 


We obtain after some easy transformations: 


t P 

'• kO 


/ _ff \ (-l)<P«k>/i /(p - |k| ) ! , (p t |l 

V ’ 2 ’ / 2<p" Ik] )/2 V(P + |k| ) ! /P- 


kl - 1) 


« I 


(38) 




, f or p - k 
even 


= 0, for p - k 

odd. 

It follows from (39) that (35) can be represented in turn as a sum of the 
terms of the form: 
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, « _ Gm' + i 




( 40 ) 


exp |i - (p - 2q) (f + oj) + (m - p + 2q)^j|* PS„,p. 2 q(c°s I) , 
q = 0, 1, 2, . . . p 

where we set 


Q = tP 

pq p 


i} 


(41) 


Using Hansen's notation: 


XT>-l.-P*2q ^ exp[i(2q - p) f]. 


we re-write (40) as: 


a; = GHI! /3'^"^^/3p+iK„ Q t; 
^•npq R 'I'mpm P‘1 '^ra 


(42) 


exp 


tnQ - 


(p - 2q) at + (m - p + 2q) 


-IV P” 

2 J P-m,p-2q 


(cos I). 


From the last expression we eliminate all terms with the period equal to 
the period of revolution of the satellite or less by performing the averaging 
over the satellite's instantaneous orbit. Making use of Hansen's form\ilas 
(Tisserand, 1889) we obtain: 


a; Q t; . 

'^>'*Pq R 'I'mpm P‘1 


(43) 


exp 


|i mf) - (p - 2q) w + (m - p + 2q)^J|* P-„, p- 2 ,(cos I) X^"^'P'^‘’(e). 
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^ere 


v-p-l, p-2q 
Ao 



exp i (p - 2q)f • dM 


(44) 


P- 1 
P - 2q| 



(1 - , 


F 



Ip- 2ql -P + 2 
^ 


|p~2q| + 1; 



The hypergeometric series in the right hand side of (44) is always a polynomial 
in e^. 

Designating the modulus of ^^mpn, ^pq ^'Cmpq argument by 

a^mpq ®uid keeping onW the real part of the disturbing function we have: 


'Vmpq 


Gm' 


j3 


I) • (43-) 


mpq 


■m , p*2 q' 


■ cos [v * T + mn - (p - 2q) o) + (m - p + 2q) - a« -i , 

2 ■tmpqj 

and die dependence on v as before is not indicated explicitly. The disturbing 
functions (43') are to be used in the actual computations of the tidal effects. Some 
particular cases of (43') deserve special attention. It follows from (44) that 


and from (43) and (45): 


v*p“l, ip 

Aq 


0 . 


Q^p = n} 

't/ffipO 'i/ttkpp 


= 0, 


(45) 
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and, consequently, 




mp 



q- I 


( 46 ) 


The terms of the type 



q 


Gm^ 

R 






-i-2 


q(COS I) • 


(47) 


cos |y»T+ mfl - - 2q) a; + (m - ^ + 2q)-^- J 

are the only terms vdiich remain in the expansion if Love coefficients are 
assumed to be constant. All other terms represent "cross-effects" be- 
tween the terms of the type (47) and deviation of Love coefficients from their 
mean values. For this reason the terms of the type (47) can be designated as 
"primary" ones. 

Among the primary terms of a special interest are the terms of the type 
^ 2 * m 2 » • . . . , + 2s). They do not contain the argument of the 

satellite's perigee) and, as a consequence, they do not change the eccentricity 
and the shape of the orbit. The same can be said about the terms ^ of 
a slightly more general type. Terms do not contain the ri^t ascension 

of the satellite's ascending node and, as a consequence, do not change the orbital 
inclination. Finally, the terms 
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( 48 ) 


^.0.2t,i R ^,0.2t.i 2« 


j(-a«“i,o cos(v ‘ T 


a™ ) 

-v.o.at.i 


V-2S-1.0 

Aq 


(1 -e^) 


2\-2»+l/2 


a- 


s, 1 - s, 1; e 


•r 


contain only the periodic terms associated with the motion of the Moon. In par- 
ticular, they carry very long period tidal effects as produced by the regression 
of the lunar node along the ecliptic. 

Setting m = q = 0 and taking 

Qoo = CC^e. 0. 0, 0) = + 1 

Into consideration, we have from (40) 


O' . CM . m' . ^io\o 

^^■Cooo " r M 2-^ + 1 


^0- 


(49) 


the terms additive to the first term GM/r In the geopotential. Theoretically, if 
not taken Into account over a long Interval of time, they will cause spurious 
small oscillations In GM, roughly one order higher than the tidal perturbations 
In the oblateness coefficient . In the motion of the satellite they can produce 
very small perturbations with the periods approximately equal to the period of 
revolution of the satellite. 
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DIFFERENTIAL EQUATIONS FOR VARIATION OF ELEMENTS 
Substituting 

fi} = n^a^BPX cos W, 

^mpq 

v^ere we set for brevity 


, (51) 

M <-mpq 

!)• (®2) 

V''»-»'»(e) (53) 

W = w *T + mO - (p - 2q) O) + (m - p + 2q)-^ - (54) 


for the disturbing function Q in the differential equations 

dsa _ 1 . ^ 

na^ vl - sin I 


(55) 


dS I _ ^ cos I 

na^ /l - sin I 



s in I 


na' 



(56) 


dSo) 

dt 


cos 1 

na* /I - e* s in I 


dn Vl - an 

na*e ‘ 


(57) 


dSe 

dt 


V I - an 

ao) ’ 


na*e 


(58) 


dSM ^ 1 - e* ^ ^ 

dt ae na aa 


(59) 
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We obtain 


dt 


= + 


1 - e* sin I 


bS?Xo cos W, 


(60) 


dSI 

dt 


n 


v^l - sin I 


[(p - 2q) cos I + m] BPX s i n W. 


(61) 


dSw _ I cos I 

” \ /rr?Tsin i ® 


X ^ ~ P ^ J B cos W 

de } 


(62) 


dSe _ n/1 - 

dt e 


(p - 2q) BPX sin W, 


(63) 


d5M r 1 - e’ dX ,, ^1 «« 

_= n — — — + 2(p t l)Xj BPcos W (64) 

Computation of generalized aasoclate Legendre functions P (cos I) and of their 

nin 

derivatives with respect to I can be accomplished using numerous recursive 
relations. The main relations are (Vilenkin, 1965): 

[(^ + n) ({ . n + 1)1''’ pL./"” 

t (n cos I - m) P^(cos I) 

s i n I 

+ [({-0X^40+ 1))'"’ I) = 0- 

I> . nco^l-” pt (cos I) (®«) 

dl ' sin I 

+ [('t - n) ('t + n ♦ 1)1*^^ 
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p|_,^(coi 1) = sin’-^ i 

(70) 

P| (cos 1) = (-1)^ • sin^ I, 

2H'. 

(71) 

I) = I)' 

(72) 

pf„(cos I) = C-ir*" pi.,-„(«s I). 

(73) 

pi„(cos I) = pJ.(cosI), 

(74) 


(75) 

Pi„(. CO. I)= (-l/-P. .„(cos I). 

(76) 
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CONCLUSION 


Balmino analysis (1974) points out the existence of the influence of the 
lateral inhomogeneities in the stricture of the Earth cm the tidal perturbations 
in the motion of a satellite. The same conclusion follows from the investiga- 
tions by Douglas et al (1974). In our opinion, in addition to analysis of pertur- 
bations of the satellite, a future planning should include the tidal observations 
on the Earth surface. With them it would be easier to obtain a proper value for 
the exterior tidal potential. The present theory points out the existence of sev- 
eral tidal long period and *'croas effects’* in the coefficients in the expansion of 
geopotential and in the motion of satellites. How long can we continue without 
including these effects and what are in fact the ’’average" elastic parameters 
which are being presently used to represent the observations of satellites ? 
These questions constitute topics for a future research. 
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